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(57) Abstract: The invention relates to a method of controlling a sound field reproduction unit (2) comprising numerous reproduc- 
tion elements (3^, using a plurality of sound information input signals (SI) which are each associated with a general pre-determined 
reproduction direction which is defined in relation to a given point (5). The invention is characterised in that it consists in: determin- 
ing parameters which are representative of the position of the elements (3 n ) in the three spatial dimensions; determining matching 
filters (A) from said spatial characteristics and said general pre-determined reproduction directions; determining control signals by 
applying the aforementioned filters to the sound information input signals (SI); and delivering control signals for application to the 
above-mentioned reproduction elements (3n). 
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(57) Abrege : Proc£d6 de pilotage d'un ensemble de restitution (2) d'un champ acoustique comportant une plurality d'elements de 
restitution (3„) a partir d'une plurality de signaux d'entree d'informations acoustiques (SI) associes chacun a une direction generate de 
restitution pre^&errninee d£finie par rapport a un point donne (5), caracterisg en ce qu'il comporte : - la determination de parametres 
representatifs de la position des Elements (3„) dans les trois dimensions de l'espace ; - la determination de filtres d'adaptation (A) a 
partir de ces caracteristiques spatiales et desdites directions generales de restitution pre<i6terminee ; - la determination de signaux de 
pilotage par l'application desdits filtres aux signaux d'entree d'informations acoustiques (SI) ; et - la delivrance de signaux de pilotage 
en vue d'une application auxdits elements de restitution (3J. 
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Procede et dispositif de pilotage d'un ensemble 
de restitution a partir d un signal multicanal. 

La presente invention concerne un procede et un dispositif de pilotage 
d'un ensemble de restitution d'un champ acoustique comportant une pluralite 
d'elements de restitution, a partir d'une pluralite de signaux acoustiques ou au- 
diophoniques assoc'ies chacun a une direction generate de restitution predeter- 
5 minee, definie par rapport a un point de I'espace donne. 

Un tel ensemble de signaux est couramment designe par I'expression 
« signal multicanal » et correspond a une pluralite de signaux, appeles canaux, 
transmis en parallele ou multiplexes les uns avec les autres, chacun destine a un 
element ou un groupe d'elements de restitution, dispose dans une direction gene- 
1 0 rale predefinie par rapport a un point donne. 

Par exemple, un systeme multicanal classique et connu sous le nom 
de « 5.1 ITU-R BF 775-1 » et comporte cinq canaux destines a des elements de 
restitutions places dans cinq directions generates predetermines par rapport a 
un centre d'ecoute, definies par les angles 0°, + 30°, -30°, +1 10° et -110°. 
1 5 Une telle disposition correspond done a la disposition d'un haut-parleur 

ou un groupe de haut-parleurs devant au centre, un de chaque c6te devant a 
droite et a gauche et un de chaque cote derriere a droite et a gauche. 

Les signaux de pilotage etant associes chacun a une direction deter- 
minee, I'application de ces signaux a un ensemble de restitution dont les ele- 
20 ments ne repondent pas a la configuration spatiale predeterminee, entratne des 
deformations importantes du champ acoustique restitue. 

II existe des systemes qui integrant des moyens de retard sur les ca- 
naux, afin de compenser au moins partiellement, I'eloignement des elements de 
restitution par rapport au centre d'ecoute. Ces systemes ne permettent cepen- 
25 dant pas de prendre en compte la disposition dans I'espace de I'ensemble de 
restitution. 

II apparatt done qu'aucun procede ou systeme existant, ne permet une 
restitution de bonne qualite a partir d'un signal de type multicanal avec un en- 
semble de restitution de configuration spatiale quelconque. 
30 La presente invention a pour objet de remedier a ce probleme, en de- 

finissant un procede et un systeme de pilotage de I'ensemble de restitution dont 
la configuration spatiale est quelconque. 
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Uinvention a pour objet un proc6de de pilotage d'un ensemble de res- 
titution d'un champ acoustique comportant une pluralite d'6Iements de restitution 
associes chacun a une direction generate de restitution predeterminee definie par 
rapport a un point donne, pour obtenir un champ acoustique restitue de caracte- 
5 ristiques specifiques sensiblement independantes des caracteristiques intrinse- 
ques de restitution dudit ensemble, caracterise en ce qu'il comporte : 

- une §tape de determination de caracteristiques au moins spatiales 
dudit ensemble de restitution, permettant la determination de parametres repre- 
sentatifs pour au moins un etement dudit ensemble de restitution de sa position 

10 dans les trois dimensions de I'espace par rapport audit point donne ; 

- une etape de determination de filtres d'adaptation a partir desdites 
caracteristiques au moins spatiales dudit ensemble de restitution et desdites di- 
rections generates de restitution predeterminee associees a ladite pluralite de 
signaux d'entree d'informations acoustiques ; 

15 - une etape de determination d'au moins un signal de pilotage desdits 

elements dudit ensemble de restitution par Tapplication desdits filtres 
d'adaptation a ladite pluralite de signaux d'entree d'informations acoustiques; et 

- une etape de delivrance dudit au moins un signal de pilotage en vue 
d'une application auxdits elements de restitution. 

20 Suivant d'autres caracteristiques : 

- ladite etape de determination de caracteristiques au moins spatiales 
dudit ensemble de restitution, comporte une sous-etape de saisie permettant de 
determiner tout ou partie des caracteristiques dudit ensemble de restitution ; 

- ladite etape de determination de caracteristiques au moins spatiales 
25 dudit ensemble de restitution, comporte une etape de calibrage permettant de 

delivrer tout ou partie des caracteristiques dudit ensemble de restitution ; 

- ladite sous-etape de calibrage comporte pour au moins un des ele- 
ments de restitution : 

- une sous-etape d'emission d'un signal specifique vers ledit au 
30 moins un element dudit ensemble de restitution ; 

- une sous-etape d'acquisition de Tonde sonore emise en re- 
ponse par ledit au moins un element ; 

- une sous-6tape de transformation desdits signaux acquis en un 
nombre fini desdits coefficients representatifs de Tonde sonore emise ; et 
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- une sous-etape de determination de parametres spatiaux et/ou 
acoustiques dudit element a partir desdits coefficients representatifs de I'onde 
sonore emise ; 

- ladite sous-etape de calibrage comporte en outre une sous-etape de 
5 determination de la position dans au moins Tune des trois dimensions de I'espace 

dudit au moins un element dudit ensemble de restitution ; 

- ladite etape de calibrage comporte une sous-etape de determination 
de la reponse en frequence dudit au moins un element dudit ensemble de restitu- 
tion ; 

10 - ladite etape de determination de filtres d'adaptation comprend : 

- une sous-etape de determination d'une matrice de decodage 
representative de filtres permettant la compensation des alterations de restitution 
dues aux caracteristiques spatiales dudit ensemble de restitution ; 

- une sous-etape de determination d'une matrice de rayonnement 
15 multicanal ideale representative des directions generates predetermines asso- 

cies a chaque signal d'informations de la pluralite des signaux d'entree ; et 

- une sous-etape de determination d'une matrice representative 
desdits filtres d'adaptation a partir de ladite matrice de decodage et de ladite ma- 
trice de rayonnement multicanal ; 

20 - ladite etape de determination de filtres d'adaptation comporte une 

pluralite de sous-etapes de calcul permettant de delivrer un ordre limite de preci- 
sion spatiale des filtres d'adaptation, une matrice correspondant a une fenetre 
spatiale representative de la repartition dans I'espace de la precision souhaitee 
lors de la reconstruction du champ acoustique et une matrice representative du 

25 rayonnement de I'ensemble de restitution, ladite sous-etape de calcul de la ma- 
trice de decodage etant realisee a partir des resultats de ces sous-etapes de cal- 
cul ; 

- les matrices de decodage, de rayonnement multicanal ideal et 
d'adaptation sont independantes de la frequence, I'etape de determination d'au 

30 moins un signal de pilotage desdits elements dudit ensemble de restitution par 
Papplication desdits filtres d'adaptation correspondant a de simples combinaisons 
lineaires suivies de retard. 

- ladite §tape de determination de caracteristiques dudit ensemble de 
restitution permet la determination de caracteristiques acoustiques dudit ensem- 
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ble de restitution et Iedit precede comporte une etape de determination de filtres 
de compensation de ces caracteristiques acoustiques, ladite etape de determina- 
tion d'au moins un signal de pilotage comprenant alors une sous-etape 
d'application desdits filtres de compensation acoustique ; 
5 - ladite etape de determination de caracteristiques acoustiques est 

adaptee pour delivrer des parametres representatifs pour au moins un element 
de sa reponse en frequence ; 

- ladite etape de determination d'au moins un signal de pilotage com- 
porte une sous-etape d'ajustement de gain et d'application de retards afin 

10 d'aligner temporellement le front d'onde des elements de restitution en fonction 
de leur distance par rapport audit point donne. 

L'invention a egalement pour objet un programme d'ordinateur com- 
prenant des instructions de code de programme pour I'execution des etapes du 
procede lorsque Iedit programme est execute par un ordinateur. 

15 L'invention a encore egalement pour objet un support amovible du 

type comprenant au moins un processeur de traitement et un element de me- 
moire non volatile, caracterise en ce que ladite m§moire comprend un pro- 
gramme comprenant des instructions de code pour I'execution des etapes du 
procede, lorsque Iedit processeur execute Iedit programme. 

20 L'invention a encore pour objet un dispositif de pilotage d'un ensemble 

de restitution d'un champ acoustique comportant une pluralite d'elements de res- 
titution, comportant des moyens d'entree d'une pluralite de signaux d'entree 
d'informations acoustiques associes chacun a une direction generate de restitu- 
tion predeterminee definie par rapport a un point donne, caracterise en ce qu'il 

25 comporte en outre : 

- des moyens de determination de caracteristiques au moins spatiales 
dudit ensemble de restitution, permettant la determination de parametres repre- 
sentatifs pour au moins un element dudit ensemble de restitution de sa position 
dans les trois dimensions de I'espace par rapport audit point donne ; 

30 - des moyens de determination de filtres d'adaptation a partir desdites 

caracteristiques au moins spatiales dudit ensemble de restitution et des direc- 
tions generates de restitution predeterminee associees a ladite pluralite des si- 
gnaux d'entree d'informations acoustiques ; et 
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- des moyens de determination d'au moins un signal de pilotage des- 
dits 6l6ments dudit ensemble de restitution par ['application desdits filtres 
d'adaptation a ladite pluralite de signaux d'entree d'informations acoustiques. 

Suivant d'autres caracteristiques de ce dispositif : 
5 - lesdits moyens de determination des caracteristiques au moins spa- 

tiales dudit ensemble de restitution comportent des moyens de saisie directe 
desdites caracteristiques ; 

- il est adapte pour etre associe a des moyens de calibrage permettant 
la determination des caracteristiques au moins spatiales dudit ensemble de 

10 restitution ; 

- lesdits moyens de calibrage comprennent des moyens d'acquisition 
d'une onde sonore comportant quatre capteurs de pression disposes selon une 
forme generate de tetraedre ; 

- lesdits moyens de determination de caracteristiques sont adaptes 
15 pour la determination de caracteristiques acoustiques d'au moins un desdits ele- 
ments de restitution dudit ensemble de restitution, ledit dispositif comportant des 
moyens de determination de filtres de compensation acoustique a partir desdites 
caracteristiques acoustiques et lesdits moyens de determination d'au moins un 
signal de pilotage etant adaptes pour I'application desdits filtres de compensa- 

20 tion acoustique ; 

- lesdits moyens de determination des caracteristiques acoustiques 
sont adaptes pour la determination de la reponse en frequence desdits elements 
de I'ensemble de restitution. 

L'invention a egalement pour objet un appareil de traitement de don- 
25 nees audio et video comportant des moyens de determination d'une pluralite de 
signaux d'entr§e d'informations acoustiques associes chacun a une direction ge- 
nerate de restitution predeterminee definie par un point donne, caracterise en ce 
qu'il comporte en outre un dispositif de pilotage d'un ensemble de restitution ; 

- lesdits moyens de determination d'une pluralite de signaux d'entree 
30 sont formes d'une unite de lecture et de decodage des disques audio et/ou video 

numeriques. 

Uinvention sera mieux comprise a la lecture de la description qui va 
suivre, donnee uniquement a titre d'exemple et faite en se referant aux dessins 
annexes, sur lesquels : 
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- la Fig.1 est une representation d'un repere spherique ; 

- la Fig.2 est un schema d'un systeme de restitution selon I'invention ; 

- la Fig.3 est un organigramme du procede de I'invention ; 

- la Fig .4 est un schema de moyens de calibrage mis en ceuvre dans 
5 le procede de I'invention ; 

- la Fig. 5 est un organigramme detaille de I'etape de calibrage ; 

- la Fig.6 est une representation simplifiee d'un capteur utilise pour la 
mise en oeuvre de I'etape de calibrage ; 

- la Fig.7 est un organigramme detaille de I'etape de determination de 
10 filtres d'adaptation ; et 

- les Figs. 8 et 9 sont des schemas de moyens de determination de si- 
gnaux de pilotage ; et 

- la Fig.10 est un schema d'un mode de realisation d'un dispositif met- 
tant en ceuvre le procede de I'invention. 

15 Sur la figure 1, on a represents un repere spherique classique, de 

maniere a preciser le systeme de coordonnees auquel il est fait reference 
dans le texte. 

Ce repere est un repere orthonormal, d'origine O et comportant trois 
axes (OX), (OY) et (OZ). 

20 Dans ce repere, une position notee x est decrite au moyen de ses 

coordonnees spheriques (r,0,0), ou r designe la distance par rapport a I'origine O, 
01'orientation dans le plan vertical et ^ I'orientation dans le plan horizontal. 

Dans un tel repere, un champ acoustique est connu si Ton definit en 
tout point a chaque instant t la pression acoustique notee p(r,0,0,f), dont la trans- 
25 formee de Fourier temporelle est notee P(r, 9, <f>,f) ou /designe la frequence. 

L'invention est fondee sur I'utilisation d'une famille de fonctions spatio- 
temporelles permettant de decrire les caracteristiques de tout champ acoustique. 

Dans le mode de realisation decrit, ces fonctions sont des fonctions di- 
tes de Fourier-Bessel spheriques de premiere espece appelees par la suite fonc- 
30 tions de Fourier-Bessel. 

Dans une zone vide de sources sonores et vide d'obstacles, les fonc- 
tions de Fourier-Bessel sont solutions de I'equation des ondes et constituent une 
base qui engendre tous les champs acoustiques produits par des sources sono- 
res situees a I'exterieur de cette zone. 



WO 2004/068463 A ^ 

T/FR2004/000115 

7 

Tout champ acoustique tridimensionnel s'exprime done par une corn- 
binaison lineaire des fonctions de Fourier-Bessel, selon Pexpression de la trans- 
formee de Fourier-Bessel inverse qui s'exprime : 

5 Dans cette Equation, les termes Pi^if) sont, par definition, les coeffi- 

cients de Fourier-Bessel du champ p(r,0,0,t), k=^- , c est la celerite du son dans 
I'air (340 ms" 1 ), ji(kr) est la fonction de Bessel spherique de premiere espdee 
d'ordre / definie par j,(x)=^ J Mn (x) oil J v (x) est la fonction de Bessel de pre- 
miere espece d'ordre v, et ;>>/"( 6>0 est I'harmonique spherique reelle d'ordre / et 
10 de terme m, avec m allant de -/ a /, definie par : 



-4= P} m \cos 6) cos( m 0) pour m >0 

-4==P,°(cos(9) pour m=0 
■y2.7U 

P} m \cos 6) sin( m <f>) pour m <0 



^4n 

Dans cette equation, les PF(x) sont les fonctions de Legendre asso- 
ciees definies par : 



15 avec Pfe) les polynomes de Legendre, definis par : 

Les coefficients de Fourier-Bessel s'expriment aussi dans le domaine 
temporel par les coefficients pi, m (t) correspondant a la transformee de Fourier 
temporelle inverse des coefficients P/, TO (/). 

20 En variante, le procede de Tinvention utilise des bases de fonctions 

s'exprimant comme des combinaisons lineaires, eventuellement infinies, de fonc- 
tions de Fourier-Bessel. 

Sur la figure 2, on a represents schematiquement un systeme de resti- 
tution dans lequel on met en ceuvre le proc6de de I'invention. 

25 Ce systeme comporte un decodeur ou adaptateur 1 pilotant un en- 

semble de restitution 2 qui comporte une pluralite d'elements 3i a 3 N , tels que 
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des haut-parleurs, des enceintes acoustiques ou toute autre source sonore ou 
groupe de sources sonores, agences de maniere quelconque dans un lieu 
d'ecoute 4. On place arbitrairement, dans le lieu d'ecoute 4, Torigine O du repere 
que Ton appelle centre 5 de I'ensemble de restitution. 
5 L'ensemble des caracteristiques spatiales, acoustiques et electrody- 

namiques est considere comme les caracteristiques intrinseques de I'ensemble 
de restitution 2. 

L'adaptateur 1 report en entree un signal SI de type multicanal compor- 
tant des informations acoustiques a restituer et un signal de definition SL compor- 
10 tant des informations representatives de caracteristiques au moins spatiales de 
I'ensemble de restitution 2 et notamment permettant la determination de parame- 
tres representatifs pour au moins un element 3 n de I'ensemble de restitution 2 de 
sa position dans les trois dimensions de Pespace par rapport au point donne 5. 

A Tissue du traitement correspondant au precede de Tinvention, 
15 l'adaptateur 1 emet & Tattention de chacun des elements ou groupes d'elements 
3i a 3a/ de I'ensemble de restitution 2, un signal sc-i a sc/vde pilotage specifique. 

Sur la figure 3, on a represents schematiquement les etapes principa- 
les du procede selon Tinvention mis en oeuvre avec un systeme de restitution tel 
que celui decrit en reference a la figure 2. 
20 Ce proced§ comporte une etape 10 de determination de parametres 

de fonctionnement, adaptee pour permettre au moins la determination des carac- 
teristiques spatiales de I'ensemble de restitution 2. 

L'etape 10 comporte une etape 20 de saisie des parametres et/ou une 
etape 30 de calibrage permettant de determiner et/ou mesurer des caracteristi- 
25 ques de I'ensemble de restitution 2. 

Dans le mode de realisation decrit, l'etape 10 comporte egalement une 
§tape 40 de determination de parametres de description des directions generates 
predetermines assoctees aux differents canaux du signal d'entree multicanal SL 
A Tissue de Tetape 10, des informations relatives au moins aux diffe- 
30 rentes directions generates predeterminees assoctees a chacun des canaux 
d'entree ainsi qu'S la position dans les trois dimensions de Tespace de chacun 
des elements ou groupes d'elements 3 n de I'ensemble de restitution 2, sont de- 
terminees. 
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Ces informations sont utilisees lors d'une §tape 50 de determination 
des filtres d'adaptation permettant de prendre en compte les caracteristiques 
spatiales de I'ensemble de restitution 2 afin de definir des filtres d'adaptation du 
signal d'entree multicanal a la configuration spatiale specifique de I'ensemble de 
5 restitution 2. 

Avantageusement, I'etape 10 permet egalement de determiner des ca- 
racteristiques acoustiques pour tout ou partie des elements 3i a 3 N de I'ensemble 
de restitution 2. 

Dans ce cas, le procede comporte une etape 60 de determination de 
10 filtres de compensation acoustique permettant de compenser I'infiuence des ca- 
racteristiques acoustiques specifiques des elements 3i a 3a/. 

Les filtres definis lors des etapes 50 et avantageusement 60, peuvent 
ainsi etre memorises, de sorte que les etapes 10, 50 et 60 ne doivent etre repe- 
tees qu'en cas de modification de la configuration spatiale de I'ensemble de resti- 
15 tution 2 et/ou de la nature du signal d'entree multicanal. 

Le procede comporte ensuite une etape 70 de determination des si- 
gnaux de pilotage sci a sca/ destines aux elements de I'ensemble de restitution 2, 
comportant une sous-etape 80 d'application des filtres d'adaptation determines 
lors de Tetape 50 aux differents canaux a(f) a c Q (t) formant le signal multicanal 
20 d'entree SI et avantageusement, une sous-etape 90 d'application des filtres de 
compensation acoustique determines lors de Tetape 60. 

Les signaux sci a scn ainsi delivres, sont appliques aux elements 3i a 
3n de Tensemble de restitution 2, afin de restituer le champ acoustique represen- 
ts par le signal multicanal d'entree SI avec une adaptation optimum aux caracte- 
25 ristiques spatiales et avantageusement acoustiques, de I'ensemble de restitution 
2. 

II apparaTt done que grace a la mise en oeuvre du procede de 
I'invention, les caracteristiques du champ acoustique restitue sont sensiblement 
independantes des caracteristiques intrinseques de restitution de I'ensemble de 
30 restitution 2 et notamment de sa configuration spatiale. 

On va maintenant decrire plus en detail les etapes principals du 
procede de I'invention. 
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Lors de I'etape 20 de saisie de parametres un operateur ou un sys- 
teme de memoire adapte peut specifier tout ou partie des parametres de calcul et 
notamment : 

- des parametres x n exprimes dans le repere sphSrique au moyen des 
5 coordonnees r ni 9* et fa , et representatifs de la position des elements 3„ par rap- 
port au centre d'ecoute 5 ; et/ou 

- des parametres H n (f), representatifs de la reponse en frequence des 
elements 3 n . 

Cette etape 20 est mise en oeuvre au moyen d'une interface de type 
10 classique telle qu'un micro-ordinateur ou tout autre moyen approprie. 

On va maintenant decrire plus en detail I'etape 30 de calibrage ainsi 
que des moyens de mise en ceuvre de cette etape. 

Sur la figure 4 on a represents le detail de moyens de calibrage. Us 
comportent un module 91 de decomposition, un module 92 de determination de 
15 reponse impulsionnelle et un module 93 de determination de parametres de cali- 
brage. 

Les moyens de calibrage sont adaptes pour etre connectes a un dis- 
positif d'acquisition sonore 100 tel qu'un microphone ou tout autre dispositif adap- 
te, et pour etre connectes tour a tour a chaque element 3 n de I'ensemble de resti- 
20 tution 2 afin de prelever des informations sur cet element. 

Sur la figure 5, on a represents le detail d'un mode de realisation de 
I'etape 30 de calibrage mise en oeuvre par les moyens de calibrage decrits pre- 
cedemment et permettant de mesurer des caracteristiques de Tensemble de res- 
titution 2. 

25 Lors d'une sous-etape 32, les moyens de calibrage emettent un signal 

specifique u n (t) tel qu'une sequence pseudo-aleatoire MLS (Maximum Length 
Sequence) a Tattention d'un element 3„. Le dispositif d'acquisition 100 regoit, lors 
d'une sous-etape 34, 1'onde sonore 6mise par Telement 3 n en reponse a la recep- 
tion du signal u n (t) et transmet / signaux cpi(t) a cpfe) representatifs de Tonde re- 

30 Que au module de decomposition 91 . 

Lors d'une sous-etape 36, le module de decomposition 91 decompose 
les signaux captes par le dispositif d'acquisition 100 en un nombre fini de coeffi- 
cients de Fourier-Bessel qi, m (t). 
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Par exemple, !e dispositif d'acquisition 100 est constitue de 4 capteurs 
de pression localises aux 4 sommets d'un tetraedre de rayon R ainsl que cela est 
represents en r6ference a la figure 6. Les signaux des 4 capteurs de pression 
sont alors notes cpi(f) a cp 4 (r). Les coefficients q 0fi (f) a g u (t) representatifs du 
5 champ acoustique capte se deduisent des signaux cp x (t) a cp 4 (t) selbn les rela- 
tions suivantes : 

g 0)0 (/) = J_ + Ci> 2(/) + CP 3 {f) + CP 4 (f) 

QlAf) = 8?? tojRf (°' Pl(/) - CP2 ^ " CP * W + CP *(/)) 

Dans ces relations CPitf) a CP 4 (J) sont les transformees de Fourier de 
1 0 cpi(0 a cp 4 (f) et 0 o> o(/) a Q x>x if) sont les transformees de Fourier de ? 0 ,o(0 a 

Lorsque ces coefficients sont definis par le module 91 , ils sont adres- 
ses au module 92 de determination de reponse. 

Lors d'une sous-etape 38, le module 92 de determination de reponse 
determine les reponses impulsionnelles hp lm (f) qui relient les coefficients de Fou- 
15 rier-Bessel qi, m (t) et le signal emis u„(t). La methode de determination depend du 
signal specifique emis. Le mode de realisation decrit utilise une methode adaptee 
aux signaux de type MLS, comme par exemple la methode de correlation. 

La reponse impulsionnelle delivree par le module 92 de determination 
de reponse est adressee au module 93 de determination de parametres. 
20 Lors d'une sous-etape 39, le module 93 deduit des informations sur 

des elements de I'ensemble de restitution. 

Dans le mode de realisation decrit, le module 93 de determination de 
parametres determine la distance r„ entre I'element 3„ et le centre 5 a partir de sa 
r6ponse hp 0 ,o(f) et de la mesure du temps mis par le son pour se propager de 
25 I'element 3„ au dispositif d'acquisition 100, grace a des methodes d'estimation de 
retard sur la reponse hpo,o(t). 

La direction (<9„,^) de I'element 3„ est deduite par calcul du maximum 
de la transformee de Fourier spherique inverse appliquee aux reponses hp 0fi (f) a 
fcpu(0 prises a I'instant t ou hp 0 , 0 (f) presente un maximum. Avantageusement, les 
30 coordonn6es e„ et sont estimees sur plusieurs instants, choisis de preference 
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autour de I'instant ou hpo,o(t) presente un maximum. La determination finale des 
coordonnees 0„ et fa est obtenue au moyen de techniques de moyennage entre 
les differentes estimations. 

Ainsi, dans le mode de realisation decrit, le dispositif d'acquisition 100 
5 est en mesure d'encoder sans ambiguTte I'orientation d'une source dans I'espace. 

En variante, les coordonnees 0 n et fa sont estimees a partir d'autres 
reponses parmi les hpi, m (t) disponibles ou sont estimees dans le domaine fre- 
quentiel a partir des reponses HP^Jj), correspondant aux transformees de Fou- 
rier des reponses hpi >m (t). 

10 Ainsi I'etape 30 permet de determiner les parametres r„, et fa. 

Dans le mode de realisation decrit, le module 93 delivre egalement la 
fonction de transfert H n {j) de chaque element 3 n , a partir des reponses hp Um {f) 
issues du module 92 de determination de reponse. 

Une premiere solution consiste a construire la reponse hp'o fi (t) corres- 

1 5 pondant a la selection de la partie de la reponse hp 0 ,o(t) qui comporte un signal 
non nul et denue des reflexions introduites par le lieu d'ecoute 4. La reponse en 
frequence H n (f) est deduite par transformee de Fourier de la reponse hp' 0 , 0 (f) pre- 
alablement fenetree. La fenetre peut etre choisie parmi les fenetres classiques de 
lissage, comme parexemple rectangulaire, Hamming, Hanning, etBlackman. 

20 Une seconde solution plus complexe consiste a appliquer un lissage 

sur le module et avantageusement sur la phase de la reponse en frequence 
HPoflij) obtenue par transformee de Fourier de la reponse hp 0 , 0 (t). Pour chaque 
frequence/ le lissage est obtenu par convolution de la reponse HP 0 , 0 (f) par une 
fenetre centree sur/ Cette convolution correspond a un moyennage de la re- 

25 ponse HP 0 ,df) autour de la frequence / La fenetre peut etre choisie parmi les fe- 
netres classiques, comme par exemple rectangulaires, triangles et Hamming. 
Avantageusement, la largeur de la fenetre varie avec la frequence. Par exemple, 
la largeur de la fenetre peut etre proportionnelle a la frequence /a laquelle est 
appliquee le lissage. Par rapport a une fenetre fixe, une fenetre variable avec la 

30 frequence permet d'eliminer au moins partiellement reflet de salle dans les hau- 
tes frequences tout en evitant un effet de troncature de la reponse HP 0 ,o(f) dans 
les basses frequences. 

Les sous-etapes 32 a 39 sont repetees pour tous les elements 3i a Z N 
de ('ensemble de restitution 2. 
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En variante, les moyens de calibrage comportent d'autres moyens 
d'acquisition d'informations relatives aux elements 3, a 3 N , tels que des moyens 
de mesure de position au laser, des moyens de traitement du signal mettant en 
ceuvre des techniques de formation de voies ou tout autre moyen approprie. 
5 Les moyens mettant en oeuvre l'etape 30 de calibrage sont constitues 

par exemple d'une carte electronique ou d'un programme d'ordinateur ou de tout 
autre moyen approprie. 

L'etape 40 permet ainsi que cela a ete dit precedemment, de determi- 
ner des parametres decrivant le format du signal multicanal d'entree et notam- 
1 0 ment les directions predetermines generates associees a chaque canal. 

Cette etape 40 peut correspondre a une selection par un operateur 
d'un format parmi une liste de formats associes chacun a des parametres memo- 
rises, et peut egalement correspondre a une detection automatique de format 
effectuee sur le signal multicanal d'entree. Alternativement, le procede est adapte 
15 pour un unique format de signal multicanal donne. Dans encore un autre mode 
de realisation, l'etape 40 permet a un utilisateur de specifier son propre format en 
saisissant manuellement les parametres decrivant les directions associees a 
chaque canal. 

II apparait que les etapes 20, 30 et 40 formant l'etape 10 de determi- 
20 nation de parametres, permettent au moins la determination de parametres de 
positionnement dans I'espace des elements 3„ de I'ensemble de restitution 2 et 
du format du signal multicanal SI. 

Sur la figure 7, on a represents un organigramme detaille de l'etape 50 
de determination des filtres d'adaptation. 
25 Cette etape comporte une pluralite de sous-etapes de calcul et de de- 

termination de matrices representatives des parametres determines prealable- 
ment. 

Ainsi, lors d'une sous-6tape 51, un parametre L, appele ordre limite 
representatif de la precision spatiale souhaitee lors de I'etape 50 de determina- 
30 tion des filtres d'adaptation, est determine par exemple de la maniere suivante : 

- le plus petit angle a min forme par une paire d'elements de I'ensemble 
de restitution 2 est calcule automatiquement au moyen d'une relation trigonome- 
trique, telle que par exemple : 

= ac °s(sin0 nl sm0„ 2 cos($,,-$ l2 ) + cos0„,cos0„ 2 ) 
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ctmin = min(a„;,„2) 

parmi I'ensemble des couples (wi, n2) tels que nl* n2; et 
- ensuite, I'ordre maximum L est determine automatiquement comme 
etant le plus grand entier respectant la relation suivant : 

L<nl a min . 

L'etape 50 de determination de filtres d'adaptation comprend ensuite 
une sous-etape 52 de determination d'une matrice W de ponderation du champ 
acoustique. Cette matrice W correspond a une fenetre spatiale W(rj) representa- 
tive de la repartition dans I'espace de la precision souhaitee lors de la reconstruc- 
tion du champ. Une telle fenetre permet de specifier la taille et la forme de la 
zone ou le champ doit §tre correctement reconstruit. Par exemple, il peut s'agir 
d'une boule centree sur le centre 5 de I'ensemble de restitution. Dans le mode de 
realisation decrit, la fenetre spatiale et la matrice W sont independantes de la 
frequence. 

West une matrice diagonale de taille (X+l) 2 contenant des coefficients 
de ponderation Wi et dans laquelle chaque coefficient Wi se trouve 2/+1 fois a la 
suite sur la diagonale. La matrice Wa done la forme suivante : 

0 



w 0 0 

0^ 



.w. 



0 

0 W L 



Dans le mode de realisation decrit, les valeurs prises par les coeffi- 
cients Wi sont les valeurs d'une fonction telle qu'une fen§tre de Hamming de taille 
de 2L+\ evaluee en /, de sorte que le parametre Wi est determine pour I allant de 
0a£. 

L'etape 50 comporte ensuite une sous-etape 53 de determination 
d'une matrice M representative du rayonnement de I'ensemble de restitution no- 
tamment a partir des parametres de position x n . La matrice de rayonnement M 
permet de deduire des coefficients de Fourier-Bessel representant le champ 
acoustique qu'emet chaque element 3„ de I'ensemble de restitution en fonction 
du signal qu'il recoil 
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M est une matrice de taille (L+l) 2 sur N, constitute d'elements M hm , m 
les indices l,m designant la iigne Hl+m et n designant la colonne «. La matrice M 
a done la forme suivante : 



10 



15 



20 



25 



•^0,0,1 -^0,0,2 M ofi-N 

^1,-1,1^1,-1,2 M,,.,,^ 

M \fi,\ -^1,0,2 -A^i.o.w 

M j.1,1 ^1.1,2 M hX , N 

• • • 

M L,-L,iM L . L>i M LrL ,„ 

• • • 

_M LyL>x M L>L>2 M L ^ N j 

Dans le mode de realisation decrit, les elements M /( , n ,„ s'obtiennent a 
partir d'un modele de rayonnement en ondes planes, de sorte que : 

M l, m ,n = yF(0n,<f>n) 

La matrice M ainsi definie est representative du rayonnement de 
I'ensemble de restitution. En particulier, M est representative de la configuration 
spatiale de I'ensemble de restitution. 

Les sous-etapes 51 a 53 peuvent etre executees sequentiellement ou 
simultanement. 

L'etape 50 de determination de filtres d'adaptation comprend ensuite 
une sous-etape 54 de prise en compte de I'ensemble des parametres du systeme 
de restitution 2 determines precedemment, afin de delivrer une matrice D de de- 
codage representative de filtres dits de reconstruction. 

En effet, les elements D n xJJ) de la matrice D correspondent a des fil- 
tres de reconstruction qui, appliques aux coefficients de Fourrier-Bessel P lm (f) 
d'un champ acoustique connu, permettent de determiner les signaux de pilotage 
d'un ensemble de restitution pour reproduire ce champ acoustique. 

La matrice de d6codage D est done I'inverse de la matrice de rayon- 
nement M. 

La matrice D est obtenue a partir de la matrice M au moyen de metho- 
des d'inversion sous contraintes faisant intervenir des parametres d 'optimisation 
supplementaires. 

Dans le mode de realisation decrit, l'etape 50 est adaptee pour realiser 
une optimisation grace a la matrice de pond6ration du champ acoustique W qui 
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permet notamment de reduire la distorsion spatiale dans le champ acoustique 
reproduit. 

Cette matrice D est delivree notamment a partir de la matrice M, selon 
I'expression suivante : 

D = (M T WM)- 1 M J W 
dans laquelle M T est la matrice transposee conjuguee de M. 

Dans le mode de realisation decrit, les matrices Met Wsont indepen- 
dantes de la frequence, de sorte que la matrice D est egalement independante 
de la frequence. Elle est constitute d'elements notes D nXm organises de la ma- 
niere suivante : 

A.o.0 A.1.-1 A.i.o A,i,r ■ ■ E>i,l,-l~ • A,i,o- • ■ A,l.£ 
A.0,0 A t i.-i A,i,o A,i,r • • A,£,-£ • • • A^.o- • • A.l,z 

: : : : : i : 

jAwA\r,i,-iAv,i,oA/,,.i' • -D NLr . L - • -D NtL>0 - • -D N>L>L _ 

L'etape 54 permet ainsi de delivrer la matrice D representative de fil- 
tres dits de reconstruction et permettant la reconstruction d'un champ acoustique 
a partir d'une configuration quelconque de I'ensemble de restitution. Grace a 
cette matrice, le precede de invention permet de prendre en compte la configu- 
ration de I'ensemble de restitution 2 et notamment de compenser les alterations 
du champ acoustique dues a sa configuration spatiale specifique. 

En variante, les parametres relatifs a I'ensemble de restitution 2 peu- 
vent etre variables en fonction de la frequence. 

Par exemple, dans un tel mode de realisation, chaque element D nJ/n (j) 
de la matrice D peut etre determine en associant a chacun des Nsignaux de pilo- 
tage une fonction de directivite D n (9,<f>j) specifiant a chaque frequence / 
I'amplitude, et avantageusement la phase souhaitee sur le signal de pilotage sc„ 
dans le cas d'une onde plane dans la direction 

On entend par fonction de directivite D n (0,<f>j) une fonction qui associe 
une valeur reelle ou complexe, eventuellement fonction de la frequence ou d'une 
plage de frequences, a chaque direction de I'espace. 

Dans le mode de realisation decrit, les fonctions de directivites sont in- 
dependantes de la frequence et notes D„(0,0). 

Ces fonctions de directivite D„(0,0) peuvent etre determinees en speci- 
fiant que certaines grandeurs physiques entre un champ ideal et le meme champ 
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reproduit par I'ensemble de restitution respectent des lois predetermines. Par 
exemple, ces grandeurs peuvent §tre la pression au centre et I'orientation du vec- 
teur Vitesse. Dans certains cas, on souhaite que seulement 3 signaux de pilotage 
soient actifs pour reproduire une onde plane. Les signaux de pilotage actifs, no- 
5 tes sc n1 a sc nZ , sont ceux qui alimentent les elements de restitution dont les direc- 
tions sont les plus proches de la direction ($0 de I'onde plane. Les elements de 
restitution actifs, notes 3m a 3„ 3 , ferment un triangle contenant la direction (0,0 
de I'onde plane. Dans ce cas, les valeurs des directivites D nl (0,<fr a D n3 {0,$) asso- 
cies aux 3 elements actifs 3„i a 3 n3 sont donnes par : 

r-V 

10 ol = WF*~t 

avec 

(sin d ni cos <j> ni sin 0„ a 
cos 4*1^ sin 0 na cos <j> na \ ( sin 9 cos <p\ 

sin0„ lS in<£ ni sin0„ 2 sin&„ sin ^ sin ^ 1 r = I sin6>sin^ I 1 = 
cos0 ni cos0„ a cos0 na J \ cos6» / 

Dans cette relation, a correspond au vecteur 
[D nl (0,0)...D n3 (0,0)] et les directions (0 n iM, (0 n2 ,fa) et (0 n3 ,fa 3 ) correspondent 
1 5 respectivement aux directions des elements 3m , 3 n2 et 3 n3 . 

On considere que les valeurs des directivites D n (0,0) correspondant 
aux elements de restitution non actifs sont nulles. 

La relation precedente est repetee pour K directions (0 k ,<f>k) d'ondes 
planes differentes. Ainsi, chacune des fonctions de directivite £>„(0,#) est fournie 
20 sous la forme d'une liste de K echantillons. Chaque echantillon est fournit sous la 
forme d'un couple D n (0 k ,fa))} ou est la direction de I'echantillon k 

et ou D„(0 k ,fa) est la valeur de la fonction de directivite associee au signal de pilo- 
tage sc„ pour la direction 

Pour chaque frequence / les coefficients D^ >m (f) de chaque fonction 
25 de directivite sont deduits des echantillons D n (0 k ,fay)}. Ces coefficients 

sont obtenus par inversion du processus d'echantillonnage angulaire qui permet 
de deduire les echantillons de la liste D n (0 k> fa))} a partir d'une fonction 

de directivite fournie sous forme de coefficients d'harmoniques spheriques. Cette 
inversion peut prendre differentes formes afin de contrdler I'interpolation entre les 
30 echantillons. 

Dans d'autres modes de realisation, les fonctions de directivite sont di- 
rectement fournies sous la forme de coefficients £>„,/,«(/) de type Fourrier- Bessel. 
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Les coefficients D„ Jim (f) ainsi determines sont utilises pour former la 

matrice D. 

L'etape 50 comporte ensuite une etape 55 de determination d'une ma- 
trice de rayonnement multicanale ideale S representative des directions genera- 
5 les predetermines associees a chaque canal du signal multicanal d'entree SI. 

La matrice S est representative du rayonnement d'un ensemble de res- 
titution ideal, c'est a dire respectant parfaitement les directions generates prede- 
terminees du format multicanal. Chaque element S t , m<q (f) de la matrice S permet 
de deduire les coefficients de Fourier-Bessel P\,Jf) du champ acoustique ideale- 
1 0 ment restitue par chaque canal c q (t). 

La matrice S est determinee en associant a chaque canal d'entree c q (f) 
et avantageusement pour chaque frequence / une figure de directivite represen- 
tative d'une distribution de sources supposees emettre le signal du canal c q (t). 

La distribution de sources est donnee sous la forme de coefficients 
15 d'harmoniques spheriques S,^, q (f). Les coefficients S^njf) sont ranges dans la 
matrice S de taille (X+l) 2 sur Q, oix Q est le nombre de canaux. 

Dans le mode de realisation decrit, l'etape de mise en forme associe a 
chaque canal c q (t) une source d'onde plane orientee dans la direction (0 q ,</> q ) cor- 
respondent a la direction (0$,$) associee au canal c q (f) dans le format multica- 
20 nal d'entree. Les coefficients Si^V) sont alors independants de la frequence, lis 
sont notes et s'obtiennent par la relation : 

Dans d'autres modes de realisation, la matrice de rayonnement ideale 
S associe une distribution discrete de sources d'ondes planes a certains canaux 
25 pour simuler I'effet d'une ceinture de haut-parleurs. Dans ce cas, les coefficients 
Sijn >q sont obtenus par sommation des contributions de chacune des sources ele- 
mentaires. 

Dans encore d'autres modes de realisation, la matrice de rayonnement 
ideale S associe certains canaux c q (t) a une distribution continue de sources 
30 d'ondes planes decrite par une fonction de directivite S q (0,<p). Dans ce cas, les 
coefficients de la matrice S sont obtenus directement par transformee de 
Fourier Spherique de la fonction de directivite S q (0,0). Dans ces modes de realisa- 
tion, la matrice S est independante de la frequence. 
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Dans d'autres modes de realisation plus complexes la matrice S asso- 
cie a certains canaux, une distribution de sources produisant un champ diffus. 
Dans ce cas, la matrice S varie avec la frequence. Ces modes de realisation sont 
adaptes aux formats multicanaux qui considerent differemment les canaux avants 
et arrieres. Par exemple, dans les applications destinees a la restitution dans des 
salles cinema, les canaux arrieres sont souvent destines a recreer une ambiance 
diffuse. 

Dans d'autres modes de realisation, la matrice S associe a certains 
canaux des sources sonores dont la reponse n'est pas plate. Par exemple, dans 
le cas ou le format multicanal associe au canal c g (t) une source d'onde plane de 
reponse en frequence H®(fl, les S&,Jf) varient avec la frequence et s'obtiennent 
par la relation : 

Si le format multicanal associe a certains canaux une superposition 
des types de distributions de sources pre-citees, les coefficients Si Jluq (f) de la ma- 
trice de rayonnement s'obtiennent par sommation des coefficients associes a 
chaque type de distribution de source. 

Enfin, I'etape 50 comporte une sous-etape 56 de determination d'une 
matrice d'adaptation spatiale ,4 correspondant aux filtres d'adaptation a appliquer 
au signal d'entree multicanal pour obtenir une restitution optimum prenant en 
compte la configuration spatiale de I'ensemble de restitution 2. 

La matrice d'adaptation spatiale A est obtenue a partir des matrices de 
mise en forme S et de decodage D au moyen de la relation : 

A =DS 

La matrice d'adaptation A permet de generer des signaux saj(f) a sadt) 
adaptes a la configuration spatiale de I'ensemble de restitution a partir des ca- 
naux cj(i) a c Q (f). Chaque element A,^(f) est un filtre specifiant la contribution du 
canal c g (f) au signal adapte sa„(t). Grace a la matrice d'adaptation A, le procede 
de I'invention permet la restitution optimum du champ acoustique decrit par le 
signal multicanal par un ensemble de restitution de configuration spatiale quel- 
conque. 

Dans le mode de realisation decrit, les matrices D et S sont indepen- 
dantes de la frequence et la matrice A egalement. Dans ce cas, les elements de 
la matrices sont des constantes notees A n>q et chacun des signaux adaptes sa](t) 
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a sattf) s'obtient par de simples combinaisons lineaires des canaux d'entree a(t) 
a c e (0, le cas 6cheant suivies de retard ainsi qu'il sera decrit ci-apres. 

Les filtres representes par la matrice A peuvent §tre mis en oeuvre 
sous differentes formes de filtres et /ou methodes de filtrage. Dans le cas ou les 
5 filtres utilises sont parametres directement avec des reponses en frequence, les 
coefficients A n>q (f) sont directement delivres par I'etape 50. Avantageusement, 
I'etape 50 de determination de filtres d'adaptation comprend une sous-etape 57 
de conversion afin de determiner les parametres des filtres pour d'autres metho- 
des de filtrage. 

1 0 Par exemple, les combinaisons de filtrage A n , q (f) sont converties en : 

- des reponses impulsionnelles finies a njq (t) calculees par 
transformee de Fourier temporelle inverse de A„ >g (f), chaque reponse im- 
pulsionnelle a n>q (f) est echantillonnee puis tronquee a une longueur propre 
a chaque reponse ; ou 
15 - des coefficients de filtres recursifs a reponses impulsionnel- 

les infinies calculees a partir des A nA {f) avec des methodes d'adaptation. 
A Tissue de I'etape 50 les parametres des filtres d'adaptation A„, g (f) 
sont fournis. 

L'etape 60 permet ainsi que cela a ete dit precedemment, de determi- 
20 ner les filtres de compensation des caracteristiques acoustiques des elements de 
I'ensemble de restitution 2 dans le cas ou des parametres relatifs a ces caracte- 
ristiques acoustiques tels que les reponses en frequence H n (f), sont determines 
lors de I'etape 10 de determination des parametres. 

La determination de tels filtres, notes HP(f) , a partir des reponses en 
25 frequence H„(J), peut etre realisee de maniere classique en appliquant des me- 
thodes d'inversion de filtres, comme par exemple I'inversion directe, les metho- 
des de decon volution, les methodes Wiener ou d'autres. 

En fonction des modes de realisation, la compensation porte unique- 
ment sur I'amplitude de la reponse ou encore sur Tamplitude et la phase. 
30 Cette etape 60 permet de determiner un filtre de compensation pour 

chaque Element 3„ de I'ensemble de restitution 2 en fonction de ses caracteristi- 
ques acoustiques specifiques. 

De m§me que precedemment, ces filtres peuvent etre mis en oeuvre 
sous differentes formes de filtres et/ou methodes de filtrage. Dans le cas oil les 
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filtres utilises sont parametres directement avec des reponses en frequence, les 
reponses HP(f) sont directement appliquees. Avantageusement, I'etape 60 de 
determination de filtres de compensation comprend une sous-etape de conver- 
sion afin de determiner les parametres des filtres pour d'autres methodes de fil- 
5 trage. 

Par exemple, les combinaisons de filtrage HfP(J) sont converties en : 

- des reponses impulsionnelles finies hfP(t) calculees par 
transformee de Fourier temporelle inverse de H$P(f) , chaque reponse im- 
pulsionnelle ffl(t) est echantillonnee puis tronquee a une longueur propre 

10 a chaque reponse ; ou 

- des coefficients de filtres recursifs a reponses impulsionnel- 
les infinies calculees a partir des H{P(f) avec des methodes d'adaptation. 

A Tissue de I'etape 60 les parametres des filtres de compensation 
HP(f) sontfournis. 

15 On va maintenant decrire plus en detail I'etape 70 de determination de 

signaux de pilotage. 

Cette etape 70 comporte une sous-etape 80 d'application des filtres 
d'adaptation representes par la matrices au signal multicanal d'entree SI corres- 
pondent au champ acoustique a restituer. Ainsi que cela a ete dit precedemment, 
20 les filtres d'adaptation A njq (f) integrent les parametres caracteristiques de 
I'ensemble de restitution 2. 

Lors de la sous-etape 80, des signaux adaptes sai(t) a sa^f) sont ob- 
tenus par I'application des filtres d'adaptation A n>q (f) aux canaux cj(t) a c Q (f) du 
signal SI. 

25 Dans le mode de realisation decrit, la matrice d'adaptation A est inde- 

pendante de la frequence et les coefficients d'adaptation A„, q sont appliques de la 
facon suivante : 

Q 

9=1 

L'adaptation se poursuit par un ajustement des gains et I'application 
30 de retards afin d'aligner temporellement les fronts d'onde des elements 3i a 3 W 
de I'ensemble de restitution 2 par rapport a I'element le plus eloigne. Les signaux 
adaptes sa\(t) a sadf) se deduisent des signaux vi(r) a v^i) selon I'expression : 
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Dans d'autres modes de realisation, la matrice d'adaptation A varie 
avec la frequence et les filtres d'adaptation A n>q (f) sont appliques de la fagon sui- 
vante : 

9=1 

avec C q (f) la transformee de Fourier temporelle du canal c q (t) et V n (f) 

defini par : 

Vn(f) = SAn W e -2*Jr»f'c 
fh 

ou &4„(/) est la transformee de Fourier temporelle de ra w (0. 
10 Selon la forme des parametres des filtres d'adaptation A n%q (f), chaque 

filtrage des canaux c q (f) par les filtres d'adaptation A^if) peut etre realise selon 
des methodes classiques de filtrage, telles que par exemple : 

- les parametres sont directement les reponses en frequence A^qif), et 
le filtrage est effectue dans le domaine frequentiel, par exemple, a I'aide des 

15 techniques usuelles de convolution par blocs ; 

- les parametres sont directement les reponses impulsionnelles finies 
tf#i,?(0i et ' e filtrage est effectue dans le domaine temporel par convolution ; ou 

- les parametres sont les coefficients de filtres recursifs a reponses 
impulsionnelles infinies, et le filtrage est effectue dans le domaine temporel au 

20 moyen des relations de recurrence. 

La sous etape 80 se termine par un ajustement des gains et 
Tapplication de retards afin d'aligner temporellement les fronts d'onde des ele- 
ments 3i a 3a/ de I'ensemble de restitution 2 par rapport a I'element le plus eloi- 
gne. Les signaux adaptes sa\(t) a sa^f) se deduisent des signaux vi(/) a v/^t) se- 

25 Ion I'expression : 

, f. max( r n ) - r n \ 
SCn(f) = fh v>,\t I 

Sur la figure 8, on a represents la structure de filtrage correspondant a 
la sous etape 80 duplication des filtres d'adaptation spatiale telle que decrite 
precedemment. 

30 Avantageusement, I'etape 70 comporte une sous-etape 90 de com- 

pensation des caracteristiques acoustiques de I'ensemble de restitution. Chaque 
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filtre de compensation HP(j) est applique au signal adapte sa n (t) correspondant 
afin d'obtenir le signal de pilotage sc n (t) de I'element 3 n , selon la relation : 

SC„(f) = SA„V) HP® 

ou SC„(f) est la transformee de Fourier temporelle de sc n (t) et ou SA„(f) est la 
5 transformee de Fourier temporelle de sa n (i). 

L'application des filtres de compensation H!P(f) des caracteristiques 
acoustiques est decrite en reference a la figure 9. 

Selon la forme des parametres de ces filtres, chaque filtrage des si- 
gnaux sa„(f) peut etre realise selon des methodes classiques de filtrage, telles 
1 0 que par exemple : 

- dans le cas ou les parametres de filtrage sont des reponses en fre- 
quence HfP(f) , le filtrage peut §tre effectue au moyen de methodes de filtrage 
dans le domaine frequentiel, comme par exemple des techniques de convolution 
par bloc ; 

15 - dans le cas ou les parametres de filtrage sont des r6ponses impul- 

sionnelles tfp(f), le filtrage peut etre effectue dans le domaine temporel par 
convolution temporelle ; 

- dans le cas ou les parametres de filtrage sont des coefficients de re- 
lations de recurrence, le filtrage peut etre realise dans le domaine temporel au 

20 moyen de filtres recursifs a reponse impulsionnelle infinie. 

Dans certains modes de realisation simplifies, le procede de invention 
ne compense pas les caracteristiques acoustiques specifiques des elements de 
I'ensemble de restitution. Dans ce cas, I'etape 60 ainsi que la sous-etape 90 ne 
sont pas realisees et les signaux adaptes scn(t) a sai^t) correspondent directe- 

25 ment aux signaux de pilotage sc, a sc N . 

Par l'application du procede de I'invention, chaque element 3i a 3 N re- 
coit done un signal de pilotage specifique sci a son et emet un champ acoustique 
qui contribue a la reconstruction optimale du champ acoustique a restituer. En 
effet, le pilotage simultane de I'ensemble des elements 3i a 3 N permet une re- 

30 construction optimale du champ acoustique correspondant au signal multicanal 
d'entree par I'ensemble de restitution 2 dont la configuration spatiale est quel- 
conque, soit encore ne correspond pas a une configuration fixe. 
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Par ailleurs, d'autres modes de realisation du procede de I'invention 
peuvent etre envisages et notamment des modes de realisation inspires de tech- 
niques decrites dans la demande de brevet en France deposee le 28 fevrier 
2002, sous le n° 02 02 585. 

Notamment, I'etape 50 de determination des filtres d'adaptation spa- 
tiale peut prendre en compte de nombreux parametres d'optimisation tel que : 

- G„{j), representatif du gabarit de I'element 3„de I'ensemble de restitu- 
tion specifiant la bande de frequence de fonctionnement de cet element ; 

- #tm»0!). representatif de la reponse spatio-temporelle de I'element 
3 n correspondant au champ acoustique produit dans le lieu d'ecoute 4 par 
I'element Z n , lorsque celui-ci recoit en entree un signal impulsionnel ; 

- W(r,f), decrivant pour chaque frequence / consideree une fenetre 
spatiale representative de la repartition dans I'espace de contraintes de recons- 
truction du champ acoustique, ces contraintes permettant de specifier la reparti- 
tion dans I'espace de I'effort de reconstruction du champ acoustique ; 

- Wtf), decrivant directement sous forme de ponderation des coeffi- 
cients de Fourier-Bessel et pour chaque frequence/consideree, une fenetre spa- 
tiale representative de la repartition dans I'espace de contraintes de reconstruc- 
tion du champ acoustique ; 

- R(j), representatif, pour chaque frequence/consideree, du rayon de 
la fenetre spatiale lorsque celle-ci est une boule ; 

- fi(f), representatif, pour chaque frequence /consideree, de la capaci- 
ty d'adaptation locale souhaitee a I'irregularite spatiale de la configuration de 
I'ensemble de restitution ; 

- {(4, m k )}(f), constituant pour chaque frequence/consideree, une liste 
de fonctions spatio-temporelles dont la reconstruction est imposee ; 

- L(f), imposant, pour chaque frequence/consideree, I'ordre limite de 
determination de filtres ; 

- RM(f), definissant, pour chaque frequence /consideree, le modele de 
rayonnement des elements 3 A a 3 w de I'ensemble de restitution 2. 

Tout ou partie de ces parametres d'optimisation peuvent intervenir lors 
de la sous-etape 54 de determination de la matrice de decodage D. Ainsi, comme 
cela est decrit dans la demande de brevet en France deposee sous le numero 02 
02 585, les parametres N lm>n {f) et RM(f) interviennent dans la sous-etape 53 de 



WO 2004/068463 A ^T/FR2004/000115 

25 

determination de la matrice de rayonnement M, les parametres Wtf), R(f) 

interviennent dans la sous-etape 52 de determination de la matrice W, les para- 
metres {(/*, m k )}(f) interviennent dans une sous-etape supplemental dans la 
determination d'une matrice F. La matrice de decodage D est alors determinee 
5 lors de la sous-etape 54, pour chaque frequence/ en fonction des matrices M, W 
et Fet des parametres G„(f) et ju(f). 

Toujours, selon la demande de brevet 02 02 585, le calcul de la ma- 
trice D peut etre effectue frequence par frequence en considerant uniquement les 
elements actifs pour chaque frequence consideree. Cette methode de determina- 

10 tion de la matrice 2> fait intervenir le parametre G n (J) et permet d'exploiter au 
mieux un ensemble de restitution dont les elements ont des bandes de frequen- 
ces de fonctionnement differentes. 

II apparaTt que la mise en oeuvre du precede de I'invention decrit ici 
est plus efficace et done plus rapide que les precedes existants et notamment 

15 que le precede decrit dans la demande de brevet francais deposee sous le nu- 
mero 02 02 585. 

En effet, pour adapter un signal multicanal comportant Q canaux a un 
ensemble de restitution comportant N elements avec une precision spatiale 
d'ordre L, il apparaTt que le precede de I'invention necessite QxN filtres 
20 d'adaptation au lieu des Q(L+1) 2 + (L+l) 2 N filtres necessaires a la mise en oeuvre 
du precede decrit dans la demande de brevet en France deposee sous le n° 02 
02 585. 

Par exemple, I'adaptation d'un signal « 5.1 ITU-R BF 775-1 » a un en- 
semble de restitution a 5 haut parleurs avec une precision d'ordre 5 necessite 25 
25 filtres au lieu de 360 filtres. 

Sur la figure 10, on a represente un schema d'un mode de realisation 
d'un appareil mettant en oeuvre le precede tel que decrit precedemment. 

Cet appareil comprend I'adaptateur 1 qui est forme d'une unite 1 10 de- 
livrant un signal multicanal telle qu'une unite de lecture de disques audio-video 
30 dite lecteur DVD 1 12. Le signal multicanal delivre par I'unite 1 10 est destine aux 
elements de I'ensemble de restitution 2. Le format de ce signal SI est reconnu 
automatiquement par I'adaptateur 1 qui est adapte pour lui faire correspondre 
des parametres decrivant la direction generate predeterminee associee a chaque 
canal du signal SI. 
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Selon I'invention, cet adaptateur 1 integre egalement une unite de cal- 
cul supplemental 1 14 ainsi que des moyens de saisie d'informations 116. 

Par exemple, les moyens de saisie 1 16 sont formes d'une interface in- 
frarouge avec une telecommande ou encore avec un ordinateur et permettent a 
un utilisateur de determiner les parametres definissant les positions dans 
I'espace des elements de restitution 3i a 3 W . 

Ces differents parametres sont utilises par le calculates 114 pour de- 
terminer la matrices definissant les filtres d'adaptation. 

Ulterieurement, le calculates 114 applique ces filtres d'adaptation au 
signal multicanal SI afin de delivrer les signaux de pilotage sc, a sc N a destination 
de I'ensemble de restitution 2. 

Bien entendu, le dispositif mettant en oeuvre I'invention peut prendre 
d'autres formes, telles qu'un logiciel mis en oeuvre sur un ordinateur ou encore 
un dispositif complet integrant des moyens de calibrage ainsi que des moyens de 
saisie et de determination des caracteristiques de I'ensemble de restitution plus 
complet. 

Ainsi, le precede peut egalement etre mis en oeuvre sous la forme 
d'un dispositif dedie a ('optimisation de systemes de restitution multicanal, exte- 
rieur a un decodeur audio-video et associe a celui-ci. Dans ce cas, le dispositif 
est adapte pour recevoir en entree un signal multicanal et delivrer en sortie des 
signaux de pilotage d'elements d'un ensemble de restitution. 

Avantageusement, le dispositif est adapte pour etre connecte au dis- 
positif 100 d'acquisition necessaire a I'etape de calibrage et/ou est muni d'une 
interface permettant de saisir des parametres, notamment, la position des ele- 
ments de I'ensemble de restitution et eventuellement le format multicanal 
d'entree. 

Un tel dispositif d'acquisition 100 peut etre relie de maniere filaire ou 
non filaire (radio, infra-rouge) et peut etre integre a un accessoire, comme une 
telecommande, ou §tre ind6pendant. 

Le procede peut etre mis en oeuvre par un dispositif integre dans un 
element d'une chame audio-video charge du traitement de signaux multicanal, 
comme par exemple un processeur ou decodeur dit « surround », un amplifica- 
teur audio-video integrant des fonctions de decodage multicanal ou encore une 
chame audio-video completement integree. 
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Le procede de I'invention peut egalement §tre mis en oeuvre dans une 
carte electronique ou dans une puce dediee. Avantageusement, il peut §tre inte- 
gre sous la forme d'un programme dans un processeur de traitement de signal 
(DSP). 

Le procede peut prendre la forme d'un programme informatique desti- 
ne a etre execute par un ordinateur. Le programme recoit en entree un signal 
multicanal et delivre les signaux de pilotage d'un ensemble de restitution even- 
tuellement integr§ a cet ordinateur. 

Par ailleurs, les moyens de calibrage peuvent etre realises en mettant 
en oeuvre un procede different de celui decrit precedemment, tel que par exem- 
ple, un procede inspire de techniques decrites dans la demande de brevet en 
France deposee le 7 mai 2002 sous le numero 02 05 741 . 
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REVINDICATIONS 

1. Precede de pilotage d'un ensemble de restitution (2) d'un champ 
acoustique comportant une pluralite d'elements de restitution (3„) a partir d'une 
pluralite de signaux d'entree d'informations acoustiques (SI) associes chacun a 
une direction generate de restitution predetermine definie par rapport a un point 
donne (5) de I'espace, pour obtenir un champ acoustique restitue de caracteristi- 
ques specifiques sensiblement independantes des caracteristiques intrinseques 
de restitution dudit ensemble (2), caracterise en ce qu'il comporte : 

- une etape (10) de determination de caracteristiques au moins spatia- 
les dudit ensemble de restitution (2), permettant la determination de parametres 
representatifs pour au moins un element (3 n ) dudit ensemble de restitution (2) de 
sa position dans les trois dimensions de I'espace par rapport audit point donne 
(5) ; 

- une etape (50) de determination de filtres d'adaptation (A) a partir 
desdites caracteristiques au moins spatiales dudit ensemble de restitution (2) et 
desdites directions generates de restitution predetermine associees a ladite plu- 
ralite de signaux d'entree d'informations acoustiques (SI) ; 

- une etape (70) de determination d'au moins un signal de pilotage 
desdits elements dudit ensemble de restitution par Implication desdits filtres 
d'adaptation a ladite pluralite de signaux d'entree d'informations acoustiques (SI); 
et 

- une etape de delivrance dudit au moins un signal de pilotage en vue 
d'une application auxdits elements de restitution (3 n ). 

2. Precede selon la revendication 1, caracterise en ce que ladite etape 
(10) de determination de caracteristiques au moins spatiales dudit ensemble de 
restitution (2), comporte une sous-etape (20) de saisie permettant de determiner 
tout ou partie des caracteristiques dudit ensemble de restitution (2). 

3. Procede selon I'une quelconque des revendications 1 et 2, caracte- 
rise en ce que ladite etape (10) de determination de caracteristiques au moins 
spatiales dudit ensemble de restitution (2), comporte une etape (30) de calibrage 
permettant de delivrer tout ou partie des caracteristiques dudit ensemble de 
restitution (2). 



10 
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4. Procede selon la revendication 3, caracterise en ce que ladite sous- 
6tape de calibrage (30) comporte pour au moins un des elements de restitution 
(3„): 

- une sous-etape Remission (32) d'un signal specifique (u n (t)) vers le- 
dit au moins un element (3 n ) dudit ensemble de restitution (2) ; 

- une sous-etape d'acquisition (34) de I'onde sonore emise en reponse 
par ledit au moins un element (3 n ) ; 

- une sous-etape de transformation (36) desdits signaux acquis en un 
nombre fini desdits coefficients representatifs de I'onde sonore emise ; et 

- une sous-etape (39) de determination de parametres spatiaux et/ou 
acoustiques dudit element (3„) a partir desdits coefficients representatifs de 
I'onde sonore emise. 

5. Precede selon I'une quelconque des revendications 3 et 4, caracte- 
rise en ce que ladite sous-etape de calibrage (30) comporte en outre une sous- 

15 etape de determination de la position dans au moins I'une des trois dimensions 
de I'espace dudit au moins un element (3 n ) dudit ensemble de restitution (2). 

6. Procede selon I'une quelconque des revendications 3 a 5, caracteri- 
se en ce que ladite etape de calibrage (30) comporte une sous-etape de determi- 
nation de la reponse en frequence (H„(f)) dudit au moins un element (3„) dudit 
ensemble de restitution (2). 

7. Procede selon I'une quelconque des revendications 1 a 6, caracteri- 
se en ce que ladite etape (50) de determination de filtres d'adaptation comprend : 

- une sous-etape (54) de determination d'une matrice de decodage (D) 
representative de filtres permettant la compensation des alterations de restitution 
dues aux caracteristiques spatiales dudit ensemble de restitution (2) ; 

- une sous-etape (55) de determination d'une matrice de rayonnement 
multicanal ideate (S) representative des directions generates predetermines as- 
socies a chaque signal d'informations de la pluralite des signaux d'entree (SI) ; et 

- une sous-etape (56) de determination d'une matrice (A) representa- 
tive desdits filtres d'adaptation a partir de ladite matrice de decodage (D) et de 
ladite matrice de rayonnement multicanal (S). 

8. Procede selon la revendication 7, caracterise en ce que ladite etape 
(50) de determination de filtres d'adaptation comporte une pluralite de sous- 
etapes de calcul (51, 52, 53) permettant de delivrer un ordre (L) limite de la preci- 
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sion spatiale des filtres d'adaptation, une matrice (W) correspondant a une fene- 
tre spatiale representative de la repartition dans I'espace de la precision souhai- 
tee lors de la reconstruction du champ acoustique et une matrice (M) representa- 
tive du rayonnement de I'ensemble de restitution (2), ladite sous-etape (54) de 
5 calcul de la matrice (2) de decodage (D) etant realisee a partir des resultats de 
ces sous-etapes de calcul. 

9. Procede selon I'une des revendications 7 ou 8, caracterise en ce 
que les matrices de decodage (D), de rayonnement multicanal ideal S et d'adap- 
tation A sont independantes de la frequence, I'etape (70) de determination d'au 

10 moins un signal de pilotage desdits elements dudit ensemble de restitution par 
I'application desdits filtres d'adaptation correspondant a de simples combinaisons 
lineaires suivies de retard. 

10. Procede selon I'une quelconque des revendications 1 a 9, caracte- 
rise en ce que ladite etape (10) de determination de caracteristiques dudit en- 

15 semble de restitution (2) permet la determination de caracteristiques acoustiques 
dudit ensemble de restitution (2) et en ce que ledit procede comporte une etape 
(60) de determination de filtres de compensation de ces caracteristiques acousti- 
ques, ladite etape (70) de determination d'au moins un signal de pilotage com- 
prenant alors une sous-etape (90) d'application desdits filtres de compensation 

20 acoustique. 

11. Procede selon la revendication 10, caracterise en ce que ladite 
etape (10) de determination de caracteristiques acoustiques est adaptee pour 
delivrer des parametres representatifs pour au moins un element (3 n ) de sa re- 
ponse en frequence (H n (f)). 

25 12- Procede selon I'une quelconque des revendications 1 a 11, carac- 

terise en ce que ladite etape (70) de determination d'au moins un signal de pilo- 
tage comporte une sous-etape d'ajustement de gain et d'application de retards 
afin d'aligner temporellement le front d'onde des elements de restitution (3 n ) en 
fonction de leur distance par rapport audit point donne (5). 

30 13 - Programme d'ordinateur comprenant des instructions de code de 

programme pour I'execution des etapes du procede selon I'une quelconque des 
revendications 1 a 12 lorsque ledit programme est execute par un ordinateur. 

14. Support amovible du type comprenant au moins un processeur de 
traitement et un element de memoire non volatile, caracterise en ce que ladite 
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memoire comprend un programme comprenant des instructions de code pour 
I'execution des etapes du procede selon I'une quelconque des revendications 1 a 
12, lorsque ledit processeur execute ledit programme. 

15. Dispositif de pilotage d'un ensemble de restitution (2) d'un champ 
5 acoustique comportant une pluralite d'elements de restitution (3 n ), comportant 

des moyens (112) d'entree d'une pluralite de signaux d'entree d'informations 
acoustiques (SI) associes chacun a une direction generate de restitution prede- 
termine definie par rapport a un point donne (5), caracterise en ce qu'il com- 
porte en outre : 

10 - des moyens (116) de determination de caracteristiques au moins 

spatiales dudit ensemble de restitution (2), permettant la determination de para- 
metres representatifs pour au moins un element (3 n ) dudit ensemble de restitu- 
tion (2) de sa position dans les trois dimensions de I'espace par rapport audit 
point donne (5) ; 

15 - des moyens (1 14) de determination de filtres d'adaptation (A) a partir 

desdites caracteristiques au moins spatiales dudit ensemble de restitution (2) et 
des directions generates de restitution predeterminee associees a ladite pluralite 
des signaux d'entree d'informations acoustiques (SI) ; et 

- des moyens (114) de determination d'au moins un signal de pilotage 

20 (sc n ) desdits elements (3 n ) dudit ensemble de restitution (2) par I'application des- 
dits filtres d'adaptation (A) a ladite pluralite de signaux d'entree d'informations 
acoustiques (SI). 

16. Dispositif selon la revendication 15, caracterise en ce que lesdits 
moyens de determination des caracteristiques au moins spatiales dudit ensemble 

25 de restitution (2) comportent des moyens (116) de saisie directe desdites carac- 
teristiques. 

17. Dispositif selon I'une quelconque des revendications 15 et 16, ca- 
racterise en ce qu'il est adapte pour etre associe a des moyens de calibrage (91, 
92, 93, 100) permettant la determination des caracteristiques au moins spatiales 

30 dudit ensemble de restitution (2). 

18. Dispositif selon la revendication 17, caracterise en ce que lesdits 
moyens de calibrage comprennent des moyens d'acquisition d'une onde sonore 
(100) comportant quatre capteurs de pression disposes selon une forme generate 
de tetraedre 
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19. Dispositif selon Tune quelconque des revendications 15 a 18, ca- 
racterise en ce que lesdits moyens de determination de caracteristiques sont 
adaptes pour la determination de caracteristiques acoustiques d'au moins un 
desdits elements (3„) dudit ensemble de restitution (2), ledit dispositif comportant 
des moyens de determination de filtres de compensation acoustique a partir des- 
dites caracteristiques acoustiques et lesdits moyens de determination d'au moins 
un signal de pilotage etant adaptes pour I'application desdits filtres de compensa- 
tion acoustique. 

20. Dispositif selon la revendication 19, caracterise en ce que lesdits 
moyens de determination des caracteristiques acoustiques sont adaptes pour la 
determination de la reponse en frequence (H n (f)) desdits elements (3 n ) de 
I'ensemble de restitution (2). 

21. Appareil de traitement de donnees audio et video comportant des 
moyens (112) de determination d'une pluralite de signaux d'entree d'informations 
acoustiques (SI) associes chacun a une direction generate de restitution prede- 
termine definie par un point donne (5), caracterise en ce qu'il comporte en outre 
un dispositif de pilotage d ? un ensemble de restitution (2) selon I'une quelconque 
des revendications 1 a 1 9. 

22. Appareil selon la revendication 21, caracterise en ce que lesdits 
moyens de d6termination d'une pluralite de signaux d'entree sont fornrtes d'une 
unite (1 12) de lecture et de decodage des disques audio et/ou video numeriques. 
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